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Understanding hadronic 7-ray emission from supernova remnants 
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Princeton University, 4 Ivy I jn > 08540 Princeton NJ, US 

We aim to test the plausibility of a theoretical framework in which the 7— ray emission detected 
from supernova remnants is of hadronic origin, i.e., due to the decay of neutral pions produced in 
nuclear collisions involving relativistic nuclei. In particular, we investigate how the nature of the 
circumstellar medium affects the evolution of a remnant and of its 7— ray emission, stressing the 
role of magnetic field amplification in the prediction of expected particle spectra. A phenomeno- 
logical scenario consistent with both the underlying Physics and the larger and larger amount of 
observational data provided by the present generation of 7— ray experiments is finally outlined and 
critically discussed. 



I. 



RAYS FROM GALACTIC SNRS 



Supernova Remnants (SNRs) have been regarded 
for many years as prominent sites for the acceleration 
of relativistic particles up to above 10 6 GeV [see e.g.[l|, 
for a review]. However, while the evidences of electron 
acceleration are clear-cut in the shape of synchrotron 
emission, a radiative evidence of the presence of non- 
thermal hadrons is much more elusive and has to be 
searched in the 7— rays resulting from the decay of 
neutral pions produced in nuclear collisions between 
cosmic rays (CRs) and the background plasma [sec 
e.g.0. 

Unfortunately, the first observations of 7— ray 
bright SNRs have demonstrated themselves to be in- 
conclusive for assessing the origin of the detected TeV 
emission, since also inverse Compton scattering of rel- 
ativistic electrons on background photons may pro- 
vide a signal in the same band. A paradigmatic case 
in this respect is RX J1713. 7-3946, the first SNR to 
be detected in 7— rays ||, whose properties have been 
accounted for in both a hadronic or a leptonic scenario 
[see e.g.[5|, and references therein]. Nevertheless, with 
the advent of the Fermi telescope, the detection of 
RX J1713. 7-3946 also in the GeV band [6] eventually 
confirmed the leptonic nature of such an emission. 

On the other hand, as showed by Q , the first strong 
radiative evidence of hadronic emission from a Galac- 
tic SNR comes from another object very recently de- 
tected by VERITAS @ and Fermi-LAT @. The rea- 
son why with broadband 7— ray observations it is pos- 
sible to infer the origin of the emission relies on the 
fact that, while hadronic 7— rays show an energy spec- 
trum parallel to the one of the parent nuclei, the 
inverse- Compton emission is instead flatter by a 0.5 
in slope. When looking at the TeV emission only, 
though, the cut-off in the electron distribution may 
show up and mimic the hadronic emission produced 
by a steep proton spectrum. Nevertheless, this effect 
cannot extend to lower energies, too: provided that 
the primary spectrum goes like ~ E~ 2 as predicted 
by first-order Fermi mechanism at strong shocks, any 
spectral slope in the GeV range larger than 1.5 may 
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FIG. 1: Spectral slope inferred for 7— ray bright SNRs, as 
a function of their age. See table 1 in [3[ for proper data 
references and ijn] for the meaning of the different curves. 



be regarded as a strong hint of hadron acceleration. 

The main properties of 7— ray-bright Galactic SNRs 
are collected in table 1 of [3[ . In particular the inferred 
slope a of the photon spectrum, both in the GeV and 
TeV bands, is reported also in fig.[T]as a function of the 
estimated age of the SNR. Black squares correspond 
to confirmed SNRs, while red diamonds correspond to 
unidentified sources which have been proposed to be 
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associated with SNRs. Green triangles, then, indicate 
SNRs whose Fermi-LAT spectrum exhibits a cut-off 
in the 1-20 GeV range: such a feature is hardly pre- 
dictable both in leptonic and hadronic models for shell 
SNRs, and may be due to the superposition of differ- 
ent sources, e.g, pulsar- wind nebulae and/or molecu- 
lar clouds (MCs). The meaning of the curves in fig. [T] 
will be explained in SjTTJ 

The most striking evidence in fig. Q] is that most 
of the detected SNRs show a photon spectral index 
larger than 2, and typically in the range a ~ 2.1-2.5. 
As outlined above, while in the TeV band one may 
be looking at the cut-off in the particle spectrum, it is 
very unlikely to account for most of these GeV spectra 
in a leptonic scenario, since an electron spectrum as 
steep as E~ 2 - 6 — E~ 3 would be required, at odds with 
what inferred, for instance, from radio measurements. 

Even more interestingly, such steep spectra is likely 
at odds also with what diffusive shock acceleration 
(DSA) at strong shocks predicts for accelerated nuclei. 
In the test-particle limit the spectral slope depends 
only on the shock compression ratio, and is a = 2 
for large Mach number shocks. However, when the 
dynamical back-reaction of the accelerated particles is 
important, namely when the pressure in CRs becomes 
a sizable fraction of the bulk one, it is possible to 
show that the expected spectrum is no longer a power- 
law, being steep er ( flatter) than E~ 2 at low (high) 
energies [see e.g. [l^, EH for some reviews]. Since in 
this context "low" means below ~ 1-10 GeV, the non- 
linear DSA theories which have been developed in the 
last three decades need to be revised under the light 
of the unprecedented wealth of data 7— ray satellites 
and Cherenkov telescopes are providing us with. 

II. MODELS AND OBSERVATIONS 

In the last few years it has become clear that at 
SNR blast waves the magnetic field can be amplified 
by a factor of few tens with respect to the mean in- 
terstellar one (1-5/iG) as a consequence of the super- 
Alfvenic streaming of the accelerated particles [see e.g. 
[T2I ] . These large magnetic fields, beside scattering CRs 
more efficiently and having a dynamical role in the 
shock dynamics [lj|, may as well significantly affect 
the shape of the spectrum of the accelerated particles. 

CRs, in fact, only feel the compression ratio of 
the waves they scatter against: if Alfven waves have 
a finite velocity v w with respect to the background 
plasma, the actual compression CRs experience is 
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where u is the fluid velocity and subscripts 1 and 
2 refer to pre- and post-shock quantities. Up- 
stream, self-generated Alfven waves travel in direc- 
tion opposite to the CR pressure gradient (v W: i = 
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FIG. 2: Slope of the energy spectrum of accelerated parti- 
cles (black lines, left axis) and self-amplified magnetic field 
upstream of the shock (red lines, right axis) as a function 
of the CR efficiency (see eqs. [2][4|). 



— va,i = — B\ / ^ An pi) , while efficient isotropization 
downstream implies v w 2 = 0. Moreover, when CR 
acceleration is efficient and the instability saturates, 
the Alfvenic Mach number ahead of the shock is cal- 
culated to be (see eq. 42 of |14|): 
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where £ cr is the CR acceleration efficiency, expressed 
here as the pressure in accelerated particles at the 
shock over the bulk pressure. Following the simple 
algebra sketched in §2 of [3|, the effective shock com- 
pression ratio at the presence of efficient CR acceler- 
ation and magnetic field amplification reads: 
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where M s is the sonic Mach number of the shock and 
7cff = 1 i s the effective adiabatic index of the 
fluid composed by CRs+thermal plasma. 

The predicted slope for the spectrum of the accel- 
erated particles is therefore simply given by 



(4) 



In fig. [2] the upstream self-generated magnetic field 
for different Mach numbers (temperature T — 10 5 K, 
number density n — 0.1 /cm 3 , velocity as in the leg- 
end) is shown on the right axis, as a function of £ cr . 
The corresponding spectral indexes expected for ac- 
celerated particles are shown on the left axis, and are 
larger and larger with increasing £ cr : this fact com- 
pletely reverses the usual prediction that the more ef- 
ficient the CR acceleration, the flatter the spectrum 
of the accelerated particles. 
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FIG. 3: Time evolution of relevant physical quantities for 
( cr = 0.1. The solid line around 1000 yr indicates the tran- 
sition between ejecta-dominated and Sedov- Taylor stages, 
while the vertical dashed lines, from left to right, mark the 
boundaries of wind zone, hot bubble and ISM, as in the 
labels (see §3 of [3| for details). 



The simple calculation above has been worked out 
within the so called two-fluid approach, in which the 
CRs are modeled as a relativistic fluid, and has there- 
fore to be checked against a kinetic modeling of parti- 
cle acceleration in which the momentum dependence 
of the spectral slope is self-consistently accounted for. 
However, preliminary results outlined in §5.1 of 14 [ 
support the scenario in which the whole particle spec- 
trum may steepen significantly when the scattering 
centers have a non-negligible velocity with respect to 
the fluid. 



III. EVOLUTION OF THE 7- RAY EMISSION 

In order to illustrate the prediction of the phe- 
nomenological model above, we consider a SNR pro- 
duced by the explosion of a very massive star which 
underwent different stages of mass ejection, and in 
particular a slow and dense red-giant wind and a fast 
and hot Wolf-Rayet wind [see e.g. [jja, and references 
therein] . 

About 80% of the Galactic SNRs are expected to be 
due to the core-collapse of massive progenitors [l6j . 
therefore the circumstellar medium their shocks prop- 
agate into should be "polluted" by winds launched in 
the pre-SN stages. The case considered here is not 
meant to describe either a given or the most common 
object: it has to be regarded only as a benchmark 
case for investigating the effects of non-trivial envi- 
ronments on the expected hadronic 7— ray luminosity 
of SNRs, actually generalizing the pioneering work of 
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The details of the semi-analytical model (usually 



referred to as thin-shell approximation) adopted to 
follow the shock propagation are widely discussed in 
§3 of What is relevant for our purposes is sum- 
marized in fig. [3] the shock (radius R s h and veloc- 
ity V s h) propagates first in the dense red- giant wind 
(density profile p oc r~ 2 ) and, for the parameters 
chosen, enters the Sedov- Taylor stage around 1000 
yr. During the Sedov stage the most energetic CRs 
are expected to escape the system (l7| carrying away 
a sizable amount of energy E esc . This energy loss 
is consistently accounted for in calculating the blast 
wave evolution. At about 4000 yr the shock enters 
the hot and rarefied bubble excavated by the fast 
(~ 2000km/s) Wolf-Rayet wind and finally, at about 
10000 yr, the shock breaks into the unperturbed in- 
terstellar medium (ISM, n = 1/cm 3 , T = 10 4 K). 

When applying the formalism of §02 the expected 
spectral slope as a function of the SNR age is plotted 
in fig. [T] for three different CR acceleration efficiencies 
£ cr = 0.01 (almost test-particle), 0.1 and 0.2, as in the 
legend. There are two very interesting points arising 
from this kind of analysis: 

• a is a rather strong function of the CR efficiency, 
but a moderate £ cr between 0.1 and 0.2 may 
account for basically all the inferred slopes; 

• for large Mach numbers the CR spectral slope 
is expected to be a function of the acceleration 
efficiency alone, and to increase when the shock 
velocity drops because of its encounter with the 
dense ISM (after 10-20xl0 4 yr). 

Clearly these qualitative results have to be checked 
against proper non-linear, kinetic, models for DSA, 
but they indeed represent a possible way to reconcile 
our current understanding of particle acceleration in 
SNRs, magnetic field amplification due to CR-induced 
streaming instability and current 7— ray detections. 
The magnetic fields needed are in fact perfectly con- 
sistent with the one inferred by X-ray measurements. 



A. The role of the circumstellar medium 

The simple model presented here also allows us to 
track the expected 7— ray luminosity of the SNR dur- 
ing its evolution, and in particular to check if the 
(rapidly growing) statistics of bright SNRs might tell 
us something about the nature of the circumstellar 
medium they are propagating into, or even give us a 
hint of when a SNR is expected to be most luminous. 
The predicted 7— ray flux above 100 MeV and 1 TeV 
is shown in fig. 0] as a function of the SNR age. 

A natural effect of the propagation through a rar- 
efied bubble (n < 0.01/cm 3 ) is in fact a net decrease 
in the expected photon flux. Provided that the statis- 
tics might not be large enough yet, and also consid- 
ered that it is quite hard to estimate correctly the 
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FIG. 4: Time evolution of the 7-ray luminosity, in the GeV 
(upper lines) and in the TeV band (lower lines). Different 
lines correspond to different CR acceleration efficiency, as 
in the legend. Vertical lines illustrate the evolutive stages 
for the case £ cr = O.f, as in figure[3] 

age of middle- age /old SNRs, the data in fig. [T] are 
consistent with a bimodal distribution, with two pop- 
ulations separated by a gap in ages around 5-10 xlO 3 
yr. Further studies may confirm or disprove such a 
hint, but indeed the present and future generation of 
7— ray telescopes could help us in understanding if 
CR acceleration occurs mainly in standard ISM or in 
superbubbles rich in Wolf-Rayet-like stars. 

B. Caveats and further comments 

The scenario depicted above may be significantly 
affected by the presence of MCs interacting with, or 



simply close to, the SNR shock. Many of the listed 
sources are in fact known to be associated with MCs 
fljl . HH which, on one hand, provide a large amount 
of targets enhancing the 7— ray luminosity but, on the 
other hand, represent a serious challenge for acceler- 
ation theories because of the issues related to shock 
propagation and particle diffusion in partially-ionized 
media. In addition, some of these sources (W28N, 
IC443, W51C and W49B, represented by green trian- 
gles in fig. [TJ exhibit a cut-off around 1-20 GeV and 
a very steep spectrum in the TeV band. While their 
position in fig. [T] is not extremely peculiar, it is likely 
that this class of sources is different from classical shell 
SNRs, therefore requiring a dedicated description. 

The model above, while retaining some key ingre- 
dients of CR-modificd shocks, does not include any 
information about the maximum energy achieved by 
particles at a given SNR age. The expected decrease 
of this quantity with the shock velocity would even- 
tually lead the 7— ray emission to fade at late evolu- 
tive stages, unless some high-energy particles remain 
trapped in the downstream. Such an effect, however, 
is very hard to quantify and is tightly related to the 
interesting problem of how accelerated particles leave 
their sources and become CRs (20j |. 

When focusing on a given object rather than trying 
a statistical approach to many Galactic remnants, it 
is indeed possible to implement more self-consistently 
the idea that magnetic field amplification may play 
a crucial role in explaining hadronic emission from 
SNRs. A nice example has been put forward by 0], 
who clearly show that Tycho's shock channels a sizable 
fraction (~ 10%) of its kinetic energy into hadrons, 
accelerated up to at least 5 x 10 5 GeV in magnetic 
fields as large as ~ 300/iG. 
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